On the shot-noise limit of a thermal current 
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The noise power spectral density of a thermal current between two macroscopic dielectric bodies 
held at different temperatures and connected only at a quantum point contact is calculated. As- 
suming the thermal energy is carried only by phonons, we model the quantum point contact as a 
mechanical link, having a harmonic spring potential. In the weak coupling, or weak-link limit, we 
find the thermal current analog of the well-known electronic shot-noise expression. 
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I. INTRODUCTION 

Just as Ohm's law, relating the electrical current to an 
applied potential, breaks down when the quantum me- 
chanical aspects of the charge carriers becomes impor- 
tant, such as in the mesoscopic regime; Fourier's Law 
of heat conduction suffers a similar fate. Mesoscopic 
phonon systems 1 provide some of the best experimen- 
tal setups to test the quantum nature of heat transport, 
such as the quantization of thermal conductance." Al- 
though, experimental demonstration of which lagged a 
decade behind that of its electronic counterpart. 3 

Following this seminal work, nanomcchanical systems 
have since seen an increased interest, experimentally 
and theoretically, from such diverse areas as quantum 
computing 4 to promising research into detecting the 
quantum mechanical zero-point motion of a macroscopic 
object. 5,6 Similar to the quantization of electrical conduc- 
tance, where each channel of a one-dimensional conduc- 
tor can contribute a quantum of electrical conductance, 
e 2 /2irh per spin, in a one-dimensional dielectric each vi- 
brational mode carries a quantum of thermal conduc- 
tance given by irk^T/6h. Of course one requirement to 
observed this quantization is a clean system with minimal 
scattering, i.e. ballistic transport. Within the Landaucr- 
Buttikcr formalism, this amounts to setting the transmis- 
sion matrix to unity for each mode. The opposite limit 
of weak transmission or strong scattering can be equally 
interesting. For instance in a system of two conductors 
separated by a thin tunnel barrier, such as a scanning 
tunneling microscope (STM), the electrical conductance, 
associated with the tunneling current, is related to the 
product of the local density of states on each side of the 
barrier. 7 In Ref. [8] a thermal analog of an STM, i.e. a 
phonon mediated scanning thermal microscope, was pro- 
posed, where the thermal conductance associated with 
the energy current between two macroscopic dielectric 
bodies held at different temperatures and connected at a 
single quantum point contact was found to be related to 
the the local phonon density of states of each reservoir. 
Similar work has been done involving the the phonon 
dominated thermal transport through more complex con- 
nections, such as molecules. 9 ' 10,11 

Here we examine the noise of a thermal current in this 
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FIG. 1: (color online) Model of two macroscopic dielectric 
bodies held at different temperatures Tl and Tr, and joined 
by a single quantum point contact. The contact can be taken 
to be (a) several atomic bonds with a spring constant K or (b) 
a small "neck" of dialectic material of length / and diameter 
d with an effective spring constant K = (nd 2 /4l)Y , where Y 
is the Young's modulus of the material. 



limit of weak transmittance, the shot-noise limit. In the 
same way the granularity of the charge carriers, in say a 
weak tunneling current, contributes to the current noise, 
the analogous behavior for phonons should be observed in 
a thermal current. 22 Experimentally, the ability to detect 
a single phonon is an ongoing area of interest. 12 

In Ref. [8] the thermal current between two insula- 
tors weakly joined by only a mesoscopic link, modeled 
as a harmonic spring, was calculated. The actual phys- 
ical link could be a few chemical bonds or even a small 
bridge of material, see Fig. [1]. The result of Ref. [8] was 
the thermal analog of the well-known tunneling current 
formula. In the following sections we examine the intrin- 
sic noise present in such a thermal current. It is assumed 
the two bodies are only weakly coupled, to lowest order 
in the coupling, this is equivalent to the shot-noise limit 
of the electronic counterpart. 



II. MODEL AND THERMAL CURRENT 

We consider the following model, which is illustrated 
in Fig. [1]: Two macroscopic dielectric bodies, labeled 
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left (L) and right (R), act as thermal reservoirs and are 
held at fixed temperatures T L and T R . The Hamiltonian 
for each side is taken, in the harmonic approximation, as 
(ft=l) 



Hj := u In a\ n a In , I = L,R 



(1) 



where a\ n and a In are phonon creation and annihilation 
operators for the left and right side, which satisfy 



[o-In^I'n'] = S ™' S 



IP 



and 



[a In ,a rn ] = [a\ n ,a\, n ,] = 0. 



(2) 



(3) 



The quantum point contact, or weak mechanical link, is 
modeled as a harmonic potential with spring constant K 



SH ;= 1 -K[ut{x ) -<(r )] : 



(4) 



where ro is the point of contact bewteen the two reser- 
voirs and Uj is the z-component (the direction normal 
to each surface) of displacement field u/(r). This model 
of the weak link assumes the comprcssional strength of 
the link dominates over other such displacements such 
as flexorial or torsional. In principle these interactions 
could also be included, this would amount to replacing 
the spring constant K with a tensor quantity coupling to 
different components of the phonon field operator. 

The displacement field of each reservoir can be ex- 
panded in terms of phonon creation and annihilation op- 
erators as 



u /( r ) := J2\J K» f /»( r ) + 4n f /n( r )] > ( 5 ) 



where f/ n (r) are the normalized vibrational eigenfunc- 
tions, and p is the mass density. 



A. Thermal Current 

Because of energy conservation and using Heisenberg's 
equation-of-motion, a thermal-current operator can be 
defined as 



Ith :=dtH K =i[H,H R }, 



(G) 



where the full Hamiltonian H = H^ + Hr + SH. Per- 
forming the commutator gives 

Ith = — y^Rnj^Rn' - A hn * , h Rn a Rn - h Rn a Rn } , (7) 



where h In := (2pco In ) 1/2 f}„, A In := h In a In + h* In a\ n , 
and {•,•} is the anticommutator. Treating the coupling 



as the perturbation; within linear response, the thermal 
current is 8 



(7th) =2kK 2 J deeNr(r ,e)N^(r ,e)[nl(e)-nl(e)], 
o 

(8) 

where nf(e) = (exp(e/fcBT J ) — I) 1 is the Bose distribu- 
tion and 

NF(t,<j) = J2 \hm(r)\ 2 S(uj - uj nI ) (9) 



is the zz component of the local spectral density. It should 
be noted that Eq. (9) is not equal to the zz component 
of the local phonon density of states tensor given by, 21 

gi{v,w) =J2fin(r)[fiM*K"-"ni), (10) 



but Nf z (r,oj) is related to the imaginary part of the re- 
tarded phonon Green's function and is the relevant quan- 
tity of interest for the present work. For clarity the su- 
perscripts zz will be dropped from here on. Eq. (8) is 
the thermal analog of the expression for an electronic 
tunneling current, Eq. (18). 



III. CALCULATION OF THE PHONON 
SHOT-NOISE 

Here we calculate the intrinsic noise 23 associated with 
a thermal current as calculated in Sec. II A. The sym- 
metrized noise is defined as 13,14 ' 15 

SthM ~\j dte^({W th (t),« th (0)}) H , (11) 

where <5ith := h\\ — \hh) H - 2i In Ref. [16] the short time, 
or high-frequency (uj — > oo), noise of a general heat cur- 
rent was studied. Here we investigate the long time or 
low- frequency (lj — ► 0) noise in the weak transmission 
limit. 

To lowest order in the interaction the noise is simply 
SM ~\J dte^({i th (t)J th (0)}) Ha , (12) 
where H ■= H L + H R , 



Tr p-P H °f) 
(O) = — - 

\ /Ho Tre-^o ' 



and 



iHnt /\^ — iHnt 



Q(t) = e lHot Oe 



(13) 



(14) 



Using Eq. (7), along with dropping anomalous terms, the 
zero-frequency component of the noise is 25 
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S th (uj = 0) = 27TA' 2 f dee 2 N h (e)N R (e){nl(e) [l + ng(e)] + ng(e) [l + n£(e)] } (15) 



or 











y (ie e 2 coth 
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iV L (e)iV R (e)[nE(e)-nl(e)]. (16) 



It is illustrative to compare Eq. (16) to the electronic expression of the zero-frequency component of the shot-noise, 

S A (u = 0) = e(I el (eV)) coth(el//3/2), (17) 

where for a tunneling current 

(l cl (eV)) =2ne\T\ 2 J2 J du Ph (ra, co - eV)p n (ra, w) [nftu - eV) - <(cj)] . (18) 



Here \T\ 2 is the transmission probability, pi(u>) is the 
electronic local density of states, and nj (u) is the Fermi 
distribution function. 

Assuming, as in most cases of interest, the phonon 
spectral density goes as a power-law at low energies, 
Nj(u>) oc to a and letting T R — ► for simplicity, the tem- 
perature dependance of the noise is given as 



5 th (w = 0) oc T 3+ 



2m 



(19) 



A. Equilibrium noise 



In the limit T T 



T R there is no net heat current; 



nonetheless, there remains thermal fluctuations given by 
Sth(w = 0) = 2fc B T 2 G th , (20) 

where 



G 



th 



lim — — 



2irK 2 / deeN L (e)N R (e) 



(21) 

is the linear thermal conductance. This is the phonon 
analog of Nyquist- Johnson noise. In an electronic system 
the Nyquist- Johnson noise is given by 



5 e i(w = 0) = 2k B TG. 



el- 



(22) 



It should be noted that, in general Eq. (20) is a univer- 
sal relation, regardless of the model used here, and is a 
consequence of the fluctuation-dissipation theorem. 



and in the low temperature limit, F c \ := S c \/I c i = e, the 
charge of the charge carrier. This has been used to mea- 
sure the fractional charge, e.g. e/3, e/5, of the quasipar- 
ticles predicted for a quantum Hall fluid. ' 18 ' 19 ' 20 

Here we determine the Fano factor for a thermal cur- 
rent. To simplfy things let Tr — > 0, thus 



Sth 
1^ 



ii 



, dee 2 N h (e)N R (e)nl( 



J oo deeiV L (e)iV R (eK( e ) 



(23) 



Again assuming a power-law form of the phonon spectral 
density and re-scaling the integrals by letting x = e/3 
gives 



f™dxx 2+2a [e x -l]~ 
j™dxx l + 2a [e x - 1] 



-k B T := C(a)k B T. (24) 



Thus the Fano factor is not a universal quantity, as in 
the electronic case, but is independent of all material 
parameters and only depends on the energy dependance 
of the spectral density. For planer surfaces 8,21 a = 1 and 
the integrals can be done analytically giving 

F th - C(l)k B T = 360 ^ (5) fc B T » 3.83 k B T, (25) 

where £(x) is the Riemann-Zeta function. One could 
loosely interpret Eq. (25) as the average energy of the 
transmitted phonons through the weak link. 



IV. DISCUSSION 



Fano Factor 



The Fano factor F, or noise-to-signal ratio, can also be 
of interest. In the case of charge shot noise, from Eq. (17) 



Besides the experimental ability to detect single 
phonons, and thus the phonon shot noise, further condi- 
tions are needed to be in the shot noise regime. Within 
the model consider here, the temperature must remain 



4 



well below any resonant modes of the weak link, also the 
link should remain in the mesoscopic regime, i.e. smaller 
than the phonon coherence length, which in itself de- 
pends on temperature. This would suggest an upper 
bound on temperatures of roughly 10 K. 

Of course phonon noise is not only of interest for the 
work presented here, but could also be used to study 
other behavior, such as demonstrating phonon bunching 
in a phonon Hanbury-Brown and Twiss experiment. 26 
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